Transition from subbarrier to deep subbarrier regimes in heavy-ion fusion reactions 
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We analyze the recent experimental data of heavy-ion fusion cross sections available up to deep 
subbarrier energies in order to discuss the threshold incident energy for a deep subbarrier fusion 
hindrance phenomenon. To this end, we employ a one-dimensional potential model with a Woods- 
Saxon internuclear potential. Fitting the experimental data in two different energy regions with 
different Woods-Saxon potentials, we define the threshold energy as an intersect of the two fusion 
excitation functions. We show that the threshold energies so extracted are in good agreement with 
the empirical systematics as well as with the values of the Krappe-Nix-Sierk (KNS) potential at the 
touching point. We also discuss the asymptotic energy shift of fusion cross sections with respect to 
the potential model calculations, and show that it decreases with decreasing energies in the deep 
subbarrier region although it takes a constant value at subbarrier energies. 

PACS numbers: 25.70.Jj,24.10.Eq 
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Heavy-ion fusion reactions at low incident energies are 
intimately related to the quantum tunneling phenomena 
of many-body systems. Because of a strong cancellation 
between the repulsive Coulomb interaction and an attrac- 
tive short range nuclear interaction between the colliding 
nuclei, a potential barrier, referred to as a Coulomb bar- 
rier, is formed, which has to be surmounted in order for 
fusion to take place. In heavy-ion reactions, because of 
a strong absorption inside the Coulomb barrier, it has 
been usually assumed that a compound nucleus is auto- 
matically formed once the Coulomb barrier has been over- 
come. The simplest approach to heavy-ion fusion reac- 
tions based on this idea, that is, a one-dimensional poten- 
tial model has been successful in reproducing experimen- 
tal fusion cross sections at energies above the Coulomb 
barrier pj. A one dimensional potential model fitted to 
reproduce fusion cross sections above the Coulomb bar- 
rier, however, have been found to underestimate fusion 
cross sections at lower energies. It has been well recog- 
nized by now that the sub-barrier fusion enhancement is 
caused by couplings of the relative motion between the 
colliding nuclei with other degrees of freedom, such as col- 
lective vibrational and rotational motions in the colliding 
nucleidi]. 

The behavior of fusion cross sections at extremely low 
energies is a critical issue for estimating reaction rates of 
astrophysical interests. One of the currents interests in 
heavy-ion fusion reactions is a steep fall-off phenomenon 
of fusion cross sections at deep subbarrier energies. Re- 
cently, fusion cross sections for several colliding systems 
have been measured down to extremely low cross sec- 
tions, up to several nb4-8]. These experimental data 
have shown that fusion cross sections fall off much more 
steeply at deep subbarrier energies as decreasing energies, 
compared to the expectation from the energy dependence 



of cross sections at subbarrier energies. Although a few 
theoretical models have been proposed d, [l(|, the origin 
for the deep subbarrier fusion hindrance has not yet been 
fully understood. 

In Refs.p. lEl Hlj]. the deep subbarrier fusion hindrance 
has been analyzed using the astrophysical S-factor. It 
has been claimed |4j,|5j,lll| that the deep subbarrier fusion 
hindrance sets in at the energy at which the astrophysical 
S-factor takes the maximum. The authors of Refs.[3, [H, 
ITTj even parametrized the threshold energy as 
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Notice that the S-factor representation provides a useful 
tool only when the penetration of the Coulomb repulsive 
potential is a dominant contribution, such as in fusion 
reactions of light systems at low energies. In fact, the 
relation between the threshold for the deep subbarrier 
hindrance and the maximum of the S-factor is not clear 
physically, and thus it is not trivial how to justify the- 
oretically the identification of the threshold energy with 
the astrophysical S-factor. Nevertheless, it has turned 
out that the threshold energy so obtained closely follows 
the values of phenomenological internucleus potentials, 
such as the Krappe-Nix-Sierk (KNS)(ll|. the BassJlJ], 
the proximity (l4j. and the Akviiz- Winther fl5j potentials, 
at the touching configuration [l6j]. This clearly implies 
that the dynamics which takes place after the colliding 
nuclei touch with each other is responsible for the deep 
subbarrier fusion hindrance, making at the same time the 
astrophysical S-factor decrease as the incident energy de- 
creases. 

In this paper, we investigate the threshold energy for 
deep subbarrier fusion hindrance using an alternative 
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method, which is physically more transparent than the 
definition with the maximum of S-factor. That is, we de- 
termine the threshold energies by fitting the experimental 
fusion cross sections in subbarrier and deep subbarrier en- 
ergy regions separately using single-channel barrier pen- 
etration model calculations, and compare them with the 
systematics given by Eq. (TTJ as well as with the touch- 
ing energy evaluated with the KNS potential. We also 
discuss the energy dependence of fusion cross sections at 
deep subbarrier energies in terms of an asymptotic energy 
shift proposed by Aguiar et «L[l"7l]. 
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FIG. 1: (Color online) Fusion excitation functions for the 
64 Ni+ 64 Ni (the upper panel) and 16 O+ 208 Pb (the lower 
panel). The solid and the dashed lines are results of single- 
channel potential model calculations which fit the experimen- 
tal data in the subbarrier and the deep subbarier energy re- 
gions, respectively. The experimental data are taken from 
Refs. Hi]. 

In order to illustrate our procedure, the upper and 
the lower panels of Fig. [1] show fusion cross sections for 
64 Ni+ 64 M and i6 O+ 208p b systemSi respectively. We first 



define the subbarrier energy region as the one in which 
fusion cross sections take between 10 -2 mb and 10° mb. 
We fit the experimental data in this energy region with 
a potential model with a Woods-Saxon potential treat- 
ing the three parameters of the potential, that is, the 
depth Vo, the radius Rq, and the surface diffuseness a, 
as adjustable parameters. To this end, we numerically 
solve the Schrodinger equation without resorting to the 
parabolic approximation^]. The fusion cross sections 
calculated in this way are shown by the solid lines in the 
figure. Of course, these calculations do not account for 
the fusion cross sections at higher energies as the chan- 
nel coupling effects are completely ignored. However, it 
is sufficient for our purpose, as we are interested only in 
the energy dependence of fusion cross sections at subbar- 
rier energies, that is, the slope of fusion excitation func- 
tions. These calculations do not reproduce the experi- 
mental data at lower energies, either. In order to obtain 
a better fit in the lower energy region, the surface diffuse- 
ness parameter has to be increased, as has been noticed in 
Refs. @, Ojj]. We then define the deep subbarrier region 
as the one in which fusion cross sections take below 10~ 3 
mb. The dashed lines in the figure show the fusion cross 
sections obtained by fitting to the experimental data in 
this energy region. See Table I for the actual values of 
the surface diffuseness parameter. From the two curves, 
we finally define the threshold energy for deep subbarrier 
fusion hindrance as the energy at which the two fusion 
excitation functions intersect with each other. 
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FIG. 2: (Color online) The threshold energy E 3 for deep sub- 
barrier hindrance for several systems, determined with the 
two slope fit to the experimental fusion cross sections, as a 
function of the parameter Z1Z2 \J A\ A2/(A\_ + A2). The solid 
curve is the empirical function given by Eq. (TTJ), while the 
stars denote the "experimental" values defined as the maxi- 
mum energy of the astrophysical S-factorsfrD]. 
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TABLE I: The threshold energy E s for deep subbarrier fusion hindrance for several systems, obtained with the two slope fit to 
the experimental fusion cross sections. a> and a< are the diffuseness parameters in the Woods-Saxon potential used to fit the 
subbarrier and the deep subbarrier regions of fusion cross sections. £ is defined as £ = Z\Z2 \/ Ai A2/(A 1 + A2), in which Zi 
and Ai (i = 1,2) are the charge and the mass numbers of the nucleus i. Ei cxp ' 1 and Ei emp ' > are the "experimental" threshold 
energy [Ly] and the empirical energies given by Eq. JTJ), respectively. Vkns is the potential energy at the touching configuration 
[l(| estimated with the KNS potential. All the energies are shown in units of MeV, while the lengths are in units of fm. 
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^Si+ b4 Ni 


1730.05 


0.71 


0.99 


46.2 


47.3±0.9 


51.3 


43.9 


16 + 208p b 


2528.55 


0.87 


0.94 


68.8 


69.6 


66.1 


70.5 


64 Ni+ 64 Ni 


4434.97 


0.76 


0.9 


88.92 


87.3±0.9 


96.1 


89.0 


60 Ni+ 89 Y 


6537.33 


0.74 


0.815 


124.5 


123±1.2 


124.5 


125.4 


90 Zr+ 89 Y 


10435.5 


0.76 


0.87 


171.8 


171±1.7 


170.3 


175.2 


90 Zr+ 90 Zr 


10733.1 


0.56 


0.76 


176.1 


175±1.8 


173.2 


179.9 


90 Zr+ 92 Zr 


10791.9 


0.53 


0.78 


171.7 


171±1.7 


173.9 


179.1 



Figure [2] shows the thr eshold energies th us obtained 
as a function of Z 1 Z 2 y/A 1 A 2 /(A 1 + A x ). The fig- 
ure also shows the threshold energy for 28 Si+ 64 Ni|20|. 
64 Ni+ 64 Nil 16 O+ 208 Pbl, 60 Ni+ 89 Y[|, 90 Zr+ 90 Zr[l|, 
90 Zr+ 92 Zr[l|], and 90 Zr+ 89 Y [^systems. For compari- 
son, the figure also shows the empirical systematics given 
by Eq. ([lj with the solid line and the "experimental" 
data defined as the maximum energy of the S-factor 

El 

by the crosses. These values are summarized in Table HI 
together with the potential energy at the touching point 
16] estimated with the KNS potential. One can see that 
the values of the threshold energy defined in our way are 
in good agreement with those defined as the maximum 
of the astrophysical S-factor as well as with the potential 
energy at the touching configuration. 

Let us next discuss briefly the asymptotic energy shift 
for deep subbarrier fusion reactions. This quantity was 
introduced by Aguiar et al. 17] as a measure of subbar- 
rier enhancement of fusion cross sections. It was defined 
as an extra energy needed to fit the experimental fusion 
cross sections with respect to a single- channel potential 
model calculation. It has been argued that the calculated 
fusion cross sections have approximately the same expo- 
nential energy dependence as the experimental data in 
the subbarrier energy region, but are shifted in energy by 
a constant amount [17| . In connection to the deep sub- 
barrier fusion hindrance, it may be interesting to revisit 
this representation. 

In order to define the asymptotic energy shift, we first 
adjust the value of Vo and i?o in the Woods-Saxon poten- 
tial, keeping the same value for the diffuseness parameter 
a as the one which has been obtained to fit to the subbar- 
rier fusion cross sections (see a> in Table I), so that the 
experimental fusion cross sections at high energies, that 
is, those above a > 100 mb, can be approximately repro- 
duced (see Fig. 13]). We then define the asymptotic energy 
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FIG. 3: (Color online) Fusion excitation functions for the 
64 Ni+ 64 Ni system. The dashed and the solid lines are results 
of single-channel potential model calculations which fit the ex- 
perimental data in the subbarrier region and at energies above 
the Coulomb barrier, respectively. The experimental data are 
taken from Refs. 0]. 



shift as a difference between the solid line in Fig. [3] and 
the experimental data for a fixed value of fusion cross 
section. Figure @] shows the asymptotic energy shift so 
extracted for several systems as a function of correspond- 
ing fusion cross section. As one can see, the asymptotic 
energy shift is nearly constant in the range of 0.1 mb 
< cr < 1 mb, in accordance to the previous conclusion by 
Aguiar et al. [17j]. However, in the deep subbarrier region, 
the asymptotic energy shift start decreasing as the fusion 
cross sections decrease, reflecting the fact that the fusion 
cross sections have a different exponential slope from that 
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in the subbarrier region, as shown in Fig. [T] 
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FIG. 4: (Color online) The asymptotic energy shift as a func- 
tion of fusion cross section for 28 Si+ 64 Ni (the filled squares), 
i6 + 208 pb ( the fiUed triangles ) ] and s 4 Ni+ 64 Ni (the filled cir- 
cles) systems. 



In summary, we have studied the energy dependence of 
heavy-ion fusion cross sections at deep subbarrier ener- 
gies using the recent experimental data. To this end, we 
employed a one-dimensional potential model. We have 
shown that the asymptotic energy shift is almost a con- 
stant in the subbarrier region, but it decreases with de- 



creasing energies in the deep subbarrier region. This is a 
clear manifestation of the hindrance phenomenon of deep 
sub-barrier fusion. In order to see at which energy the 
deep subbarrier hindrance takes place, we estimated the 
threshold energy with a two-slope fit procedure. That is, 
we defined the threshold energy as an intersect of two 
fusion excitation functions, which fit the experimental 
fusion cross sections either in the subbarrier energy re- 
gion or in the deep subbarrier energy region. We have 
shown that the threshold energies so defined are in good 
agreement with those estimated from the maximum of 
astrophysical ^-factor. 

The definition for the threshold energy proposed in this 
paper is complementary to the one with the maximum of 
astrophysical S-factor. As we have shown in this paper, 
both the definitions provide a similar value of threshold 
energy as the potential energies at the touching config- 
uration. This strongly suggests that the dynamics after 
the touching plays an important role in deep subbarrier 
fusion reactions, changing the exponential slope of fusion 
cross sections and at the same time making the astrophys- 
ical S-factor take the maximum, although it is an open 
question why and how the dynamics after the touching 
leads to the maximum of astrophysical S-factor. 
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